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^ ' Abstract 



We compute the angular two-point correlation functions of the gamma-ray bursts at cos- 
^ ■ mological distances. Since the gamma-ray burst emission mechanism is not yet established, 
we simply assume that the gamma-ray burst sources are associated with high-redshift galaxies 
in some way. Then on the basis of several simple models for the evolution of galaxy spatial 
correlations, we calculate the amplitude of angular two-point correlation functions on scales 
appropriate for the Compton Gamma Ray Observatory data, we find that in most cases the 
predicted correlations are difficult to detect with the current data rate and the angular reso- 
lution, but models in which the bursts preferentially occur at relatively low redshift (z £ 0.5) 
predict correlation amplitudes on 9 ~ 5° which will be marginally detectable with the Gamma 
Ray Observatory data in several years. If future observations detect a signal of angular corre- 
lations, it will imply useful information on the correlation of galaxies at high redshifts provided 
that the gamma-ray bursts are cosmological. 

PACS number(s): 95.85.Pw, 98.62.Py, 98.80.Es 
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1. Introduction 



The origin of the gamma-ray bursts (GRBs) remains one of the challenging puzzles in high- 
energy astrophysics. Although several models to account for the violent energy production 
have been proposed so far, none of them seems to be successful in explaining various observed 
features of the GRBs including typical light curves, energy spectra, source distributions, in a 
consistent and convincing manneilil. Nevertheless the high degree of isotropy in the angular 
distribution of GRBs observed by the BATSE (Burst and Transient Source Experiment) on 
the Gamma Ray Observatory (GRO) can be interpreted as a strong case for the cosmological 
origin whatever the emission mechanism isS In particular, a V/V max test is a very stringent 
argument in favor of cosmological origini. 

If the GRBs are cosmological, the most natural idea is that the GRB sources are associated 
with galaxies at redshifts around and up to z = (1 ~ 2)@ i. If this is the case, the angular 
distribution of the GRBs should trace that of distant galaxies, provided that the GRB event 
rate be constant in time and insensitive to the host galaxies. Therefore the cosmological model 
for the GRBs implies that the accumulating number of the GRBs detected by GRO should 
eventually reveal the underlying angular correlation which is expected for the distant galaxies. 

In the present paper, we employ several simple models for the evolution of galaxy correla- 
tion functions, and compute the angular two-point correlation functions of the GRBs assuming 
that the GRB event rate is constant. We find that some models predict correlation ampli- 
tudes on 9 ~ 5° which will be marginally detectable with the GRO data in several years. 
The predicted amplitude would become significantly larger if we assume that the GRBs are 
preferentially associated with distant rich clusters as preliminarily reported by Cohen, Kolatte 
and Pirani. 

The rest of the paper is organized as follows. In §2 our model for the evolution of galaxy 
correlations and the GRBs is presented. §3 describes the predictions of the angular two- 
point correlation wgrb{0) for the GRBs, and examines the dependence on the several model 
parameters. Finally §4 discusses the detectability of the correlation, and the implications of 
the present results. 



2. Models 

Throughout the present paper, we adopt a simple hypothesis that GRBs are of cosmological 
origin. To be more specific, we assume that they are associated with distant galaxies and that 
their event rate is independent of host galaxies, for simplicity. With these assumptions, our 
present model is basically specified by two quantities; one is the probability of the GRB event 
per galaxy per unit time, 4>(z), and the other is the spatial two-point correlation function of 
galaxies, £, gg (x, z). 

Since no specific model for the GRB sources is available which enables to predict the form 
of <p(z), we consider two simple cases. Model A assumes that the GRB event rate is constant 
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from some critical redshift z c until the present: 



<j> A {z{t)) 



< z < z c 

Z r <Z 



On the other hand, our model B assumes that the GRBs occur for a finite period Az centered 
on z c : 

(h B (z(t)) = ' Zc ~ Az / 2 - Z - Zc + Az / 2 (2) 

V 1 K)) \ ; otherwise ' 1 ; 

where 0o A and 4>o B are constant and independent of the redshift and host galaxies. 

The other important ingredient in computing the angular correlation of the GRBs is 
£ gg (x,z). Once dark matter and cosmo logical parameters are specified, it is fairly straightfor- 
ward to theoretically predict the two-point correlation functions of dark matter. The two-point 
correlation function of galaxies, however, is quite difficult to estimate due to various compli- 
cated physical processes involved in galaxy formation. Therefore, we adopt a power-law fit to 
the present galaxy-galaxy correlation functiorM 

£ gg (x,z = 0) = (x/x )-\ x = 5.4/1"^, 7 = 1-77, (3) 

where h is the Hubble constant Hq in units of 100 km/sec/Mpc. As for the time evolution, 
we assume that ^gg (x, z) increases the amplitude independently of x: 

^g (X, Z) = g 2 (z, Qq, \ )£gg(x, Z = 0) (4) 

where and Ao denote the cosmological density parameter and a dimensionless cosmological 
constant at present, and g(z, Qq, Ao) is a growth rate of the correlation function. As a growth 
rate g(z, f2o, A ), we use either the linear growth rate D\(z, flo, Xo)/D 1 (0, Q , A ), or the power- 
law parameterization (1 + ^)~( 3 + e_ '>)/ 2 . The parameter e is introduced in order to incorporate 
the possible nonlinear growth in a very crude manner. If e = 0, the clustering measured in 
proper coordinates is not changing, and if e = —3 + 7 = —1.23, that measured in comoving 
coordinates is not changing. Linear growth rate in the Einstein - de Sitter models corresponds 
to e = 7 — 1 



3. Angular two-point correlation function of GRBs 

3.1 basic equations 

We consider the Robert son- Walker metric for the background universe: 

dx 2 



ds 2 = -dt 2 + a 2 (t) 



l-Kx 2 



+ x 2 d9 2 + x 2 sin 2 



(5) 
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where a(t) is the cosmic scale factor, K is the spatial curvature constant which is related to 
the cosmological parameters as follows: 



^2 tt2 

K= a J& 



(flo + Ao-1), 



(6) 



where do = a(t ) is the scale factor at the present time to, and c is the speed of light. The 
universe is assumed to be dominated by non-relativistic matter. In order to examine the 
effect of the background universe model, we compute the angular correlation functions of the 
GRBs, wgrb(9), both in open universes without cosmological constant (fl < 0, A = 0) and 
in spatially flat universes (f2 < 0, A = 1 — £Iq)- 

Using the small separation approximation (e.g., §56 in Ref. 5 ) with the above assumptions, 
wqrb{6) is related to £ S9 as 



W G rb{0) 



F 2 (x) 



/OO 
(lu£gg(X, Z) 
-OO 



o Fix 



-i 2 



(7) 



where u = F(x)yr 2 /a 2 — x 2 9 2 , and F(x) = yl — Kx 2 . In the above expression, the redshift 
z should be regarded as a function of the comoving coordinate x through the null geodesic 
equation. More explicitly, they are related to each other as 



for A = 0, and 



H a x _ 2[(O - 2)(1 + n z) 1/2 + 2 - fi + n z 
c ~ fig(l + z) 

dz' 



(9) 



for A = 1 — flo- 

For the power-law model for £ gg (x,z) (eq.[I|), the integration over u in eq.(^) can be 
performed and w GRB {6) in this model reduces t 



A 



a e 1 - 1 , 



OO /£' 

dx 



5-7 



F{x 



-<p 2 (z)(i + zy g 2 (z,n ,\ ) 



dx^$ 
o F(x) 



1 2 



(10) 

(11) 



where if 7 is a product of Gamma functions; H y = r(l/2)T ((7 — l)/2) /T(j/2). The cosmo- 
logical model enters through the volume factor F, the relation between x and z (eqs.J§ and 
@), and g(z,Q ,\ ). 
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3.2 the Einstein - de Sitter model 



In the Einstein - de Sitter model with the power-law £ gg (x,z), the amplitude A of w GRB (9) 
(eq.jllll) can be expressed in terms of the hypergeometric function 2 Fi. 

Hr n (y^f(y c /2) ; (model A) 

A = H.i^y-^- x y°rf(y + /2) - y^/(y-/2) , (12) 

c 6 - t [ (g _ y 3 )2 ; ( model B ) 

/(y) = 2^1(6-7,27-26-6,7-7^), (13) 
y c = 2[l-(l + ^ c )- 1 / 2 ], (14) 
y± = 2[l-(l + z c ±A^/2)- 1 / 2 ]. (15) 

When Az < z c C 1, the above expressions reduce to 

9 

; (model A) 

A ~ /-y-(^)" x ^ . (16) 

; (model B) 

Fi gure 1 plots the amplitude of w GRB at 9 = 5° function of z c in the Einstein - 
Sitter model. The angular scale 9 = 5° roughly corresponds to the angular resolution for 
the GRO. Thus mainly we refer to the amplitude at this scale in the analysis below. The 
angular correlation w(5°) decreases with z c since the number of accidental neighbors from the 
foreground and background increases with z c and thus the angular correlation is reduced. The 
difference of evolution of correlation function of galaxies appears at z c > 1. At a given z, the 
spatial correlation function for smaller e is larger than that for larger e. Thus, the angular 
correlation function decreases with e, and this effect appears at high redshift. In Figure 2, we 
show a contour of w GRB (5°) for model B on a (z c ,Az/z c ) plane. 




3.3 Open and spatially-flat cosmological models 



When the background cosmological model is general or the £, gg (x, z) is not simply assumed to 
be a power-law, one has to numerically evaluate the integral of eq.(^). The linear growth rates 
of fluctuations, D\, are explicitly given by 



h b! 3/2 

b x = (Qo 1 -l)a/a 
for Aq = models, and by 



In 



(l + fol )l/2_ 6l l/2 



In 



1.37 



2 + b 2 '\ 



dbo' 



0.2b 2 + 0.0446 2 
1 + 0.3286o L37 



b 2 = 2 1 /3(fi-i _ l)V3 a/ao 



(17) 
(18) 

(19) 
(20) 
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for Ao = 1 — models. For Ao = 1 — ^o models, we use an empirical fitting formula (|T9| ) in 
the present computation!. 

Fi gures 3a and 3b illustrate w GRB (5°) as a function of Qq for models A and B, respectively 
(with z c = 0.5, and Az/z c = 0.1). When the growth rate is given independently of the value of 
Qq (as in our power-law models for g), the angular scale 9 corresponds to the smaller physical 
length in the larger Q model. Thus the amplitude of w GRB at a given 9 becomes larger for 
larger Q since it is mainly contributed by the spatial correlation £ gg on smaller length scales. 
This is not the case when the growth rate g depends on Qq- Since we normalize the amplitude 
of £, gg at the present epoch, lower Q models, especially Ao = cases, should have higher 
amplitude of £ gg in the past (since linear density fluctuations do not grow for z ^ CIq 1 .) This 
is why the amplitude of w GRB decreases with increasing Q (for Ao = cases) when we adopt 
linear growth rate for g(z, tto, Ao)- In any case, the variation of the predicted amplitude of 
w GRB for different is very small, and at most within a factor of 2. 



3.4 CDM spectrum 



So far we adopted power-law models (eq. |§) for the spatial correlations of galaxies. As an 
example based on fairly specific theoretical models, let us consider the cold dark matter (CDM) 
models assuming that galaxies trace mass (dark matter). In this case, we adopt £(x,z) gg 
computed from the CDM power spectrum P(k, z) in linear theory: 



£,gg\%-> z) 



sin kx 



k 2 dk. 



2n 2 J y 1 ' kx 

As for P(k, z), we use a fitting formula for the primordial Harrison-Zel'dovich spectrumlll: 



(21) 



P(k,z) = CoD^k 



In (1 + 2.34g) 
2.34g 



1 + 3.89g + (16. lg) 2 + (5.46g) 3 + (6.71?)' 



-1/2 



(22) 



where q = fc/(f2o/i 2 Mpc *), and the amplitude Co is fixed so that the top-hat filtered mass 
fluctuation becomes unity at r = 8/i _1 Mpc (see, e.g. Ref.8). 

We evaluate ^ gg {x, z = 0) and w GRB {9) numerically in the CDM model, which are shown in 
Figures 4 and 5 for model A in the Einstein-de Sitter model (Qq = 1.0, Ao = 0.0 and h = 0.5). 
The amplitude of w GRB (9) for the CDM model rapidly decreases if the angle is larger than the 
angle spanning the horizon scale at the equality epoch which locates at z ~ z c . This is because 
the amplitude of £ 
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x, z = 0) for the CDM model decreases rapidly outside the horizon scale 
at the equality epoch. Incidentally this may provide useful information of z c and Qq, if the 
angular scale of the zero-point of w GRB {9) is detectable^. 



4. Discussion and conclusions 



It is important to ask here whether or not the amplitude of w GRB (9) predicted in the previous 
section is detectable. This is checked as follows. The expected number of GRB pairs at 
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separation 6 to 6 + 56 is 

(DD(0))60 = ^£^£££[1 + Waw (9)]2Tr989, (23) 

2 47T 

where N GRB is the number of observed GRBs homogeneously surveyed over the 4V steradian. 
A statistically significant estimate of w GRB (8) requires that the number of pairs in excess of 
random (the second term in RHS of eq.||23||) should be at least larger than the shot noise 
fluctuation given by 

(ADD(6))56 = J Ngrb N T B 2^659. (24) 
V 2 Air 

This implies that given w GRB (8), 6, and 89, one needs 

The lower bound on N GRB for successful detection of w GRB (5°) is listed in the last row in Table 
1. Taking into account the BATSE detection rate of GRBs ~ 1 per day!, the amplitude of w GRB 
in most models is rather small for successful detection. Nevertheless it is marginally detectable 
in some models (for instance, model A with z c < 0.3, model B with z c = 1.0, Az/z c = 0.1 and 
e = —3 + 7) with GRBs observed by BATSE in several years. 

This difficulty in detecting w GRB is partly ascribed to the low angular resolution of the 
currently available 7-ray detectors. In fact, it would be much easier to detect w GRB at small 
angular separations for a given number of sources; w GRB {6) oc for power- law model (see 
Fig. 4). The correlation of GRBs at 1° is about three times as large as that at 5°, and the 
necessary number of GRBs reduces factor of 1.4. In this respect, radio, optical and X-ray 
follow-up observations, if successful, would provide us very valuable information on w GRB . 

Recently, Cohen, Kolatte and PiranH reported a possible association of GRBs with rich 
clusters. As long as neutron stars are involved in GRB events, this association would be 
difficult to explain. Therefore their report is a serious challenge for the GRB mechanism 
with neutron stars. On the other hand, if their claim turns out to be real, the amplitude of 
Wgrb would become significantly larger than computed in the previous section since clusters 
of galaxies are known to be much more strongly clustered than galaxies; depending on the 
richness of clusters which GRBs are associated with, the amplitude would be enhanced by 
even an order of magnitude. Then w GRB (5°) would be detectable in a few years for most 
evolutionary models considered here. 

In summary, we employ several simple models for the evolution of galaxy correlation func- 
tions, and compute the angular two-point correlation functions of the GRBs assuming that 
the GRB event rate is constant. On the basis of the plausible evolutional models for galaxy 
clustering up to z & 1, the predicted amplitudes of w GRB at 5° scale are generally too small 
to be detected with the current 7-ray instruments, although some (optimistic) models pre- 
dict marginally detectable values. Such small amplitudes are mainly due to our simplifying 
assumption that the GRBs faithfully reflect the clustering at high redshift galaxies. Without 
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any other convincing theory, we believe that this assumption is most plausible and realistic. 
Nevertheless the amplitudes of the angular correlation functions may become significantly 
larger than our predictions either if the GRBs are associated with highly biased species of 
galaxies (including a possible association of rich clusters) or if the galaxy clustering around 
z £ 1 is much stronger than what standard cosmo logical models predict. Obviously either 
case seems rather unrealistic, but cannot be ruled out until specific GRB model is established. 

Therefore even if no positive signature of w GRB around 5° scales is detected from future 
data in several years, this is completely consistent with their cosmological models. In turn, 
in case that w GRB are detected, the result provides us useful information on the correlation of 
galaxies at high redshifts provided that the GRBs are cosmological. 

One of authors (S.K.) acknowledges the support by the Fellowships of Japan Society for 
the Promotion of Science for Japanese Junior Scientists. This research was supported in part 
by the Grants-in-Aid by the Ministry of Education, Science and Culture of Japan (05640312, 
06233209). 
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Table 1. Summary of w GRB (5°). 



model 


Z c 


Az/z c 


^0 


Ao 


e 




N G RB 


A 


0.5 


- 


1.0 


0.0 





2.9 x 10" 3 


7900 


A 


1.0 


- 


1.0 


0.0 





9.7 x 10~ 4 


24000 


B 


0.5 


0.1 


1.0 


0.0 





1.7 x 10~ 2 


1400 


B 


1.0 


0.1 


1.0 


0.0 





6.2 x 10~ 3 


3700 


A 


0.5 


- 


1.0 


0.0 


-3+7 


4.3 x 10" 3 


5300 


A 


1.0 


- 


1.0 


0.0 


-3+7 


1.9 x 10~ 3 


12000 


B 


0.5 


0.1 


1.0 


0.0 


-3+7 


2.7 x 10~ 2 


800 


B 


1.0 


0.1 


1.0 


0.0 


-3+7 


1.5 x 10" 2 


1500 


A 


0.5 


- 


0.2 


0.0 





2.5 x 10~ 3 


9200 


A 


1.0 


- 


0.2 


0.0 





7.5 x 10~ 4 


31000 


B 


0.5 


0.1 


0.2 


0.0 





1.3 x 10" 2 


1700 


B 


1.0 


0.1 


0.2 


0.0 





4.2 x 10~ 3 


5500 


A 


0.5 


- 


0.2 


0.0 


-3+7 


3.7 x 10~ 3 


6200 


A 


1.0 


- 


0.2 


0.0 


-3+7 


1.5 x 10^ 3 


15000 


B 


0.5 


0.1 


0.2 


0.0 


-3+7 


2.2 x 10- 2 


1100 


B 


1.0 


0.1 


0.2 


0.0 


-3+7 


9.8 x 10~ 3 


2300 


A 


0.5 




0.2 


0.8 





1.9 x 10~ 3 


12000 


A 


1.0 




0.2 


0.8 





5.1 x 10~ 4 


45000 


B 


0.5 


0.1 


0.2 


0.8 





9.3 x 10~ 3 


2500 


B 


1.0 


0.1 


0.2 


0.8 





2.6 x 10~ 3 


8800 


A 


0.5 




0.2 


0.8 


-3+7 


2.9 x 10~ 3 


7900 


A 


1.0 




0.2 


0.8 


-3+7 


1.0 x 10~ 3 


23000 


B 


0.5 


0.1 


0.2 


0.8 


-3+7 


1.5 x 10~ 2 


1500 


B 


1.0 


0.1 


0.2 


0.8 


-3+7 


6.2 x 10~ 3 


3700 
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Figure captions 

Figure 1. The angular two-point correlation functions of GRBs at 6 = 5° as a function of z c 
in the Einstein - de Sitter model. The thick solid, dotted and thin solid lines represent 
results for e = 0, e = — 3 + 7 and linear growth rate, respectively, (a) model A; (b) 
model B with Az/z c = 0.1. 

Figure 2. The contourlines of wgrb{5°) in the (z c , Az/z c ) plane for model B (the Einstein 
- de Sitter model). The thick solid, dotted and thin solid lines represent results for 
e = 0, e = — 3 + 7 and linear growth rate, respectively. 

Figure 3. The angular correlation w GRB (5°) plotted against fi - Thick lines indicate results 
for models without cosmological constant, while thin lines for spatially flat models 
(Ao = 1 — fio)- Solid, dotted, and dashed lines correspond to e = and —3 + 7, and 
linear growth rate, respectively, (a) model A with z c = 0.5, (b) model B with z c = 0.5 
and Az/z c = 0.1. 

Figure 4. The two-point correlation function of galaxies plotted against separation x (the 
Einstein-de Sitter model). Solid line corresponds to the CDM model (f2 = 1-0, A = 0.0 
and h = 0.5), while dotted line corresponds to a power-law model. 

Figure 5. The angular two-point correlation function of GRBs plotted against angular 
separation 9. Thick and thin lines indicate results for CDM model and for power-law 
model. In both cases, we adopt model A for <p{z) with z c = 1.0 and 0.1, and linear 
growth rate in the Einstein-de Sitter model for g. 
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This figure "figl-2.png" is available in "png" format from: 
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